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We perform a thorough study of the theoretical consistency of recently proposed, viable, quadratic modifica- 
tions of gravity that are functions of the the Gauss-Bonnet invariant, regarding the stability of their perturbations 
around vacuum, maximally symmetric spaces of constant curvature. We pay special attention, in particular, 
to the investigation of pathological instabilities associated with the occurrence of propagating spin-0 tachyon 
modes, and with the development of a graviton ghost. The latter effect is associated with the known "Ricci 
stability" issue, well studied in /(_R)-theories. Within quadratic modifications of gravity it is discussed for the 
first time. Special attention is paid to the requirement of non-negativity of the effective gravitational coupling, 
which warrants that the graviton is not a negative-norm state. It is demonstrated that, several theories that pass 
the cosmological as well as the solar system tests, have to be rule out on the basis of the unavoidable character 
of these pathological instabilities. 
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I. INTRODUCTION 

Attempts to modify the Einstein-Hilbert (EH) action of ge- 
neral relativity (GR) have been motivated by a number of 
reasons. In particular, renormalization at one-loop demands 
that the Einstein-Hilbert action be supplemented by higher or- 
der curvature terms 111] . Besides, when quantum corrections or 
string theory are taken into account, the effective low energy 
action for pure gravity admits higher order curvature inva- 
riants \^. More recently it has been suggested that the present 
cosmic speed-up could have its origin in - among other possi- 
bilities - corrections to the GR equations of motion, generated 
by non-linear contributions of the scalar curvature R in the pu- 
re gravity Lagrangian of f{R) theories yj-ja]. Next in degree 
of complexity are the so called F{X^ Y, Z)-theories ll6l- [l2ll . 
where, for simplicity of writing, we have introduced the field 
variables X = R- the curvature scalar, Y = i^^^i?'^", and 



^auaX^ 



vg\ 



The gravitational spectrum of the lineari- 



zation of these theories consists of a massless spin-2 graviton 
plus two additional gravitational propagating degrees of free- 
dom: a massive spin-0 excitation, and a massive spin-2 propa- 
gating mode. The latter appears to be a ghost mode associated 
with the Weyl curvature invariant C^ = C^^<,aC^'""^ ^^. 
Notwithstanding, there are ways to overcome (or at least to 
smooth out) the consequences of the would be massive spin-2 
ghost mode as, for instance, in case F{X., Y, Z) ~ F{X, Q), 
where G ^ X^ - 4Y + Z is the Gauss-Bonnet (GB) invariant. 
It has been demonstrated that quadratic modifications of 
gravity where the GB invariant enters as a function f{Q) ad- 
ded to the gravitational action: F{X, Q) = X + .fJG), ™ght 
represent a candidate for effective dark-energy yl [ij, fl4ll . 
Examples of such quadratic modifications are those where po- 
wers of the GB invariant are involved. In Re. il4ll. for instan- 



ce, models of the kind /(t?) = /ol^P were considered. In that 
reference, the authors showed that, when arbitrary functions 
of the GB-invariant are added to the GR action, the resulting 
theories are endowed with a quite rich cosmological structure: 
they may naturally lead to an effective cosmological constant, 
quintessence or phantom cosmic acceleration, with a possibi- 
lity for the transition from deceleration to acceleration. It has 
been also demonstrated that these theories are viable, since 
they comply with the Solar System constraints. Specific pro- 
perties of f{G) gravity in a de Sitter universe, such as dS and 
SdS solutions, their entropy and its explicit one-loop quanti- 
zation were also studied in Ill4ll . Additionally, the issue of a 
possible solution of the hierarchy problem in modified gravi- 
ties was addressed too. A related, yet simpler, class of models 
is given by the following action: ' 



S^ 



16nG 



N 



-{x+fiGn, 



(1) 



where Gn is the Newton's constant and the parameter /i can 
be either positive or negative. For positive values of the para- 
meter ?7, the above model can be consistent with solar-system 
tests for n < 0,074, if the GB term is to be the cause of the 
present speedup of the cosmic expansion fllSll . For n < 
the model ([U is not cosmologically viable due to the existen- 
ce of separatrices between the radiation-dominated, and dark 
energy-dominated stages of the cosmic evolution 01611 . 

Recently several new F{X, ^)-models have been proposed 
lUTl UM that are peculiar functions of the GB invariant, such 
as arctan, exp, etc. In 111711 . three particular models have been 
proposed (F{X, g) = X + f{g)): 
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' This is a particular class of models f{G) 
invariant is always positive. 



fi\g\" 114], when the GB- 



f{g) = A^ arctaii (^^ - aA^^, (3) 



f{G) = A^^ln 



cosh I I 



-aA^^, 



(4) 



where A, a, and /3 are positive constants. Meanwhile, in Ref. 
lUal the following models have been explored: 



f{g) = a\/^expf - 



2/3 



long-distance modification of gravity^ develops a ghost gra- 
viton, and also, a tachyon instability associated with negati- 
ve mass squared of the spin-0 propagating perturbation ild\ ? 
Hence this theory has to be rule out as inconsistent. This de- 
monstrates the need for a careful investigation of the possible 
(2) pathological instabilities that can be associated with the linea- 
rization of the above quadratic modifications of gravity ([T]|7]) 
around vacuum, maximally symmetric spaces of constant cur- 
vature. 

Aim of the present paper is, precisely, to carefully investi- 
gate the vacuum structure of the proposed F{X, t/)-theories 
by expanding the corresponding Lagrangian densities about 
vacuum, maximally symmetric spaces of constant curvature. 
Then we look at the spectrum ofperturbations of the metric by 
comparing with known results ^^. We pay special attention 
to checking the models regarding absence of Ricci instabili- 
ties, as well as spin-0 tachyon instabilities, and, also, to chec- 

(5) king non-negativity of the effective gravitational coupling, an 
implicit requirement that is not always carefully checked. Ric- 

(6) ci instability is usually associated with the development of a 
ghost graviton and, although it has been well established for 

(7) /(X)-theories, within F{X, Y, Z) gravity theories, as long as 
we know, this issue has not been discussed before. 



where, as before, a and /3 are constants. The above models 
have been shown to be cosmologically viable and compatible 
with solar system tests. In Ref. Ill9ll . for instance, the models 
(l2]|7| were compared to combined data-sets of supernovae, 
baryon acoustic oscillations, and constraints from the CMB 
surface of last scattering. The authors found that these models 
provide data-fits that are close to those of the LCDM concor- 
dance model. It is suggested that higher order gravity models 
that pass the mentioned tests, need to be compared to cons- 
traints from large scale structure and full CMB analysis 11 1911 . 
A similar, seemingly positive situation, arises with theories 
of the kind ([T]). Actually, a generalization of the theory depic- 
ted by the action dTJ for negative /i = — /2^~"'" and n < 0, 
which is singular in the limit of vanishing curvature, has been 
discussed in details in references lilQ, 1II1I (see also 11211 ). It 
has been demonstrated, in particular, that although gravity is 
modified at large distances, one recovers an acceptable New- 
tonian limit at small distances [1 11]. It has been discussed the- 
re that, with the values of the parameter p. necessary to ex- 
plain the cosmic acceleration {p. ~ Hq), the correction to the 
Newtonian potential is negligible for solar system scales but 
becomes important at galactic lengths. In consequence, alt- 
hough having a negligible effect on the solar system, this kind 
of modifications would have implications for the dark mat- 
ter problem lllOl [Till . The above behavior can be explained 
as follows. Even if one chooses a theory ([T]i, which is spin-2, 
ghost- free (it corresponds to the choice a = c, b ^ —4c in the 
model of Ref. lIlOlUlll ). there is an additional spin-0 degree of 
freedom in the resulting theory. It then arises that the mass of 
this extra excitation depends on the background in such a way 
that it effectively decouples when close to any matter source 
(the Sun, galaxies, etc.), but in vacuum it is of the order of the 
Hubble scale. However, a careful study of the vacuum structu- 
re of the corresponding linearized theory shows that the above 



II. VACUUM INSTABILITIES 

Stability issues are central in the study of higher order mo- 
difications of general relativity, since these are plagued by se- 
veral kinds of instabihties, some of which are catastrophic, 
leading to subsequent ruling out of the corresponding theo- 
ries. Amongst these is the fundamental Ostrogradski instabi- 
lity based on the powerful no-go theorem of the same name 
1I21II : "there is a linear instability in the Hamiltonians asso- 
ciated with Lagrangians which depend upon more than one 
time derivative in such a way that the dependence cannot be 
eliminated by partial integration". This result is general and 
can be extended to higher order derivatives in general. As a 
consequence, the only Ostrogradski-stable higher order mo- 
difications of Einstein-Hilbert action are those in the form of 
an f{X) function 1I21I1 . This result alone might rule out any 
intent to consider quadratic modifications such as the ones of 
interest here. However, the subject is subtle and, in the last 
instance, consideration of such theories can shed more light 
on the stability issue. 

In order to discuss such subtle issues, usually, one expands 
the action of the theory in the neighborhood of background 
(vacuum) spaces of constant curvature. If one considers gene- 
ral actions of the kind 



^ Recall that we are considering the particular case given by (TJ with both 

fi < and n < negative quantities. 
' Here we have in mind perturbations around vacuum maximally symmetric 

spaces of constant curvature. 



d'^x 



16t:Gn 



F{X,Y,Z), 



(8) 



2A; 



/O - ^0 /g/6 + 2Xq fgg 



(13) 



then, Taylor expanding the function F{X, Y, Z) around ma- 
ximally symmetric spaces of constant curvature X = Xq, 
y = Yo = ^o/4, Z = Zo = X^/6 (the subscript "0"means 
the given invariant is evaluated at the vacuum, constant cur- 
vature background), and keeping terms up to 0(3), one gets 
(see S for details): 



F(X,Y,Z) = X + aX + ^X^ + ^C^ , (9) 
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where 



C^ = C^.^xCf^"^^ ^Z-2Y+-X^ 



and 



m^ = 



4Xo'/<" 



GQ 



^'^^afgg 



(14) 



is the mass squared of the spin-0 perturbation propagating in 
the background space. 

As seen from ( fT2l i. the linearized F{X, g) -theor y ca n be 
recast into the form of an /(A')-theory (see also ll20ll ). so 
that, at linearized level, it is also Ostrogradski-stable. In 
what follows, under the term "Ostrogradski-stable", just linear 
Ostrogradski-stability is to be meant. 

The constant 



is the Weyl curvature invariant, and the coefficients are 



a=l-4Xlf^g, 



(15) 



X^ 

X = Fq - XqFx - YqFy - ZqFz + -z-Fxx , 



Px - ^oFxx : P = "^^xx + 2-Fy + 2i^2 



l^F^Y 



4F" 



(10) 



Now, since we will be interested in theories of the kind 
F[X, r, Z) = F{X, g) ^ X + f{g), then it will be con- 
venient to make the following replacements in the above ex- 
pansion: 



^x 



XX 



l + fgQ%,F^ 
iggiSx) 






g^z^ 



fggiSxf + fgSxx^ 



(11) 



where, as before, the suffix "0" means that the given mag- 
nitude is evaluated at background curvature Xq, and we ha- 
ve taken into account the fact that, for F{X, (7)-theories of 
the form: F{X, Q) = X + f{g), Fx = 1, Fxx = 0, and 
Fg = fg. Besides, since gx = 2X, gxx == 2, C^y = -4, 
and gz = 1, then we can write the expansion of the above 
F{X, tJ)-theories around vacuum of constant curvature Xq, 
in the following form: 



in ( fT2l) plays an important role since it multiplies the gravi- 
tational action, thus modifying the gravitational coupling. In 
fact, by substituting the linearized function ( fT2b back into the 
action ^ one gets that 



Sn 



a 
2^ 



d'^x\/\g\ 



-2A + X- 



6m^ 



:X' 



so that the effective gravitational coupling is given by 



BttG, 



eff 



K 

a 



1 



4^oV<° 



(16) 



gg 



If a given F{X, tJ)-theory fails to generate a non-negative ef- 
fective gravitational coupling, then the graviton is a negati- 
ve norm state (a ghost), and the theory has to be rule out as 
unphysical. Within /(X)-theories fulfillment of this require- 
ment (non-negativity of Ge//) is known as "linear stability". 
We have to differentiate Ge// from the effective Newton's 
"constant" SttG^ = k^/Fx, which is multiplying the right- 
hand-side (RHS) of the Einstein's field equations: 



G 



fj.iy 



8^G^^^(tM+T-), 



(17) 



F{X,g) = ,fo--X^fl + 2X^f^g 

+ (1 - 4X^f")X + 2XlfggX^ + 0(3) 



which can be derived from ([8]) by varying with respect to the 

metric 16']. In dlTl i. T^w is the stress-energy tensor for matter, 
while 



The linearized function F{X, Q) can be further written in a 
form that is convenient to compare it with known results when 
the metric is perturbed around the vacuum solution ||3-S ll21 '■ 



F{X,g)^a[-2K + X 



Qm? 



X' 



(12) 



where we have defined the effective cosmological constant 



2rpcur 



-g^,{F - XFx) - (.g^.n - V^V,)F 



X 



— 2{FYR'^^j^Rcry+FzR\apiJiR "^ '' i,) — 9 fj.iy'^ a'^ \ {Fy R"^ ) 

-aiFyR^,) + 2W,Wx{FyR\^R\) + 2FzR\^,j'), 

where D = .g'^^V^Vi, is the D'Lambertian, and we have 
used the usual representation for symmetrization: T^jj^^-f = 



{T^u + Tyi^i)/2. This form of writing of the field equations 
is particularly useful for cosmological applications. The tra- 
ce of equation ( fTTj i generates an additional constraint on the 
curvature, which, for vacuum spaces (T*^™^ = 0) of constant 
curvature can be written as: 



Theories that meet the above bound are said to be Ricci sta- 
ble 1I22II24I1 . Since the mass squared of the scalar degree of 
freedom in linearized /(X) -theories is given by 



Jx 



^ofxx 



ill 



XX 



^^0 ^0^ X ^^o^'y ^-^o^z — 0. (18) ^jj(j since linear stability of the scalar perturbation requires 



As a check of consistency, notice that by taking into account 
the trace equation ( fTSl l, which when evaluated at background 
constant curvature Xq, in terms of the ^-derivatives of f{G), 
amounts to: 



^0 + 2/0 - -X^f'^ - , 



(19) 



then the effective cosmological constant computed with the 
help of ( fT3] l, yields A = Xi^/A, as it should be for a back- 
ground space of constant curvature. 

If in equations ( fTTj i one replaces Fx ^^ 1 + 2X fg, etc., 
then, one can write the effective Newton's constant in the fo- 
llowing form: 



SttG" 



N 



1 + 2A/e 



(20) 



Note that, since equations ( [TtI i hold in general, the magnitudes 
in ( l20l i have not been evaluated at constant curvature. Howe- 
ver, in order to compare with other bounds, we can indeed 
evaluate at constant X = Xq. 

Another kind of catastrophic instability, the so called "Ric- 
ci instability" 12211 . also known as Dolgov-Kawasaki instabi- 
lity 12311 (see also yl |24[]), is associated, precisely, with the 
behavior of the effective Newton's constant G^ . Although 
the issue has been well explored within /(Ar)-theories, for 
F(X, Y, Z)-theories a detailed study of this issue is lacking. 
While in f (X)-theones Ricci stability is closely related to the 
absence of spin-0 tachyon instability, for F{X, Y, Z) gravity 
theories these are separate requirements. To understand the is- 
sue one might rely on the following reasoning line. Following 
the arguments in 112211 . suppose that the magnitude of the ef- 
fective Newton's constant, which variation is governed by 



dX 



8^Px 



F- 



XX 



2K^fg+2X^fgg) 

87r(l + 2Xfg)^ 



>0 



increases as the curvature X increases. Hence, at large cur- 
vature the effect of gravity becomes stronger. Additionally, 
increasing X generates larger values of the curvature through 
the trace equation ( fTSl l. The combination of the above mecha- 
nisms conspires to destabilize the theory: a (initially) small 
curvature grows without limits and the system runs away. In 
order to avoid this kind of unwanted instability one has to re- 
quire that 



dC 



eff 

N 



dX 



<0 



fg + 2X^fgg > . 



that /^ — Xo/^jf > 1I25II . hence, absence of spin-0 tach- 
yon instability within /(X)-theories requires that f'^^ > 
(the theory has to be Ricci stable). This is why, within f{X)- 
theories Ricci stability is closely related with absence of a sca- 
lar ghost."* When dealing with F{X, Y, Z)-theories the situa- 
tion is a bit different. Actually, due to the definition of the 
mass squared of the spin-0 perturbation in Eq. ( fT4l i. the theory 
can be Ricci stable fg + 2X^fgg > 0, and, at the same ti- 
me, the scalar perturbation can be linearly stable {Geff > 0) 
1 — 4Xq /gg > 0, and yet it can develop a tachyon instability 
as long as f^g < 0. 

The instability caused by the presence of a spin-0 tachyon 
degree of freedom, can be catastrophic either since the time 
scale for this instability to manifest is estimated to be of the 
order of cic rrin , and it might develop as quickly as in a time 



Summarizing our discussion. The relevant stability requi- 
rements imposed on theoretically consistent quadratic modifi- 
cations of gravity of the kind F{X, Q) = X + f{G), are the 
following 112011 : 

■ Ostrogradski Stability: The linearized F{X, ^)-theory 
should be expressible as an equivalent /(X)-theory.^ 
Theories of the kind F{X,Q) are obviously spin-2 
ghost-free, and their linearizations can be written in the 
form of an /(X)-theory. Hence, the models of inter- 
est here ([T]|7|, are free of the Ostrogradski instability, at 
least at the linearized level. 

■ Non-negativity of the (inverse) Effective Gravitational 
Coupling: 



a 



1 - 4^oVgg 



>G 



Fulfillment of this constraint warrants that the graviton 
is not a ghost. 

Absence of Tachyon Instability: 



-4Xg/ga 
12^0 /ee 



>0 



"* Used in this context the meaning of "ghost" state is misleading, since, in 
field theory a ghost is usually associated with a state of negative "kinetic" 
energy (otherwise a state of negative norm), while a state with negative 
mass squared is usually associated with a tachyon degree of freedom. 

' Stated in this form Ostrogradski stability implies also, absence of spin-2 
Weyl ghost propagating modes. 



Once non-negativity of the effective gravitational 
coupling is verified, the requirement of absence of spin- 
tachyon instability amounts just to fgg > 0. Fulfill- 
ment of the latter constraint means that the scalar degree 
of freedom is not a destabilizing tachyon. 

■ Ricci Stability: 

/° + 2X^f^g > . 

This requirement guarantees that the graviton does not 
become a ghost. 

Note that the bound established on viable quadratic theories 
of gravity by the dynamical systems study, according to which 
stability of the de Sitter point requires that fgg > |17y, is 
consistent with the requirement of absence of tachyon insta- 
bility (third item above) once non-negativity of the effective 
gravitational coupling is achieved. As already seen the latter 
is an independent requirement. 

We want to recall that spin-2 ghost-free theories are also 
Ostrogradski stable in the sense discussed before. Notice that 
for an /(X) -theory, the last three requirements are not inde- 
pendent: Non-negativity of Gg// (linear stability) and Ricci 
stability, together imply that the scalar mode is not a tachyon. 

Now we are in position to discuss the stability of theories 
given by ([l]|7]), regarding the last three requirements above. 



III. INSTABILITIES IN THE MODEL F{X, G) = X + ^iG" 
The function F{X, Q) in the action ([Hi is 

F{x, g) = x + ^iG'' , 

where, for generality of the analysis we keep arbitrary the sign 
of the constants /i and n, i. e., /i e 5R, n e 5R. The derived 
magnitudes fg ~ nfiQ"^^, fgg ~ n{n — 1)/^CJ"^^, when 
evaluated at the background value Xq (recall that Go ~ ^o/Q) 
are given by: 



Taking into account the above relationships we get that: 



JG 



nfJ- ^2n-2 fO _6n(n-l)^ 2«-4 



6" 



-X' 



7 Jl 



QQ 



& 



-xi 



(21) 



The value of the background's constant curvature Xq is de- 
termined by means of the trace equation (fT9] l. which, in the 
present case can be written as 

X, + 2/0 - ^fg = ^ [6" + 2^(1 - n)X'o''^'] = . 

As seen, the latter requirement amounts to an algebraic equa- 
tion^ which can be solved, in principle, to get Xq as a function 
of the overall parameters: 



X' 



QT, 



2(n- 1)^ 



4^o'/< 



3^0 



72n fO ^1^ 

l^n, Jgg — ^3 , 



Jg + 2Xq jgg — 



3n(12n- 11) 
(n-l)Xo ■ 



Note that the stability bounds discussed above for the model 
([T]i, do not depend on the parameter /i. Therefore, the sign 
of the quadratic modification of GR: f{G) oc G^\ does not 
play role in the stability of its linearization around maximally 
symmetric background spaces of constant curvature. Let us 
write the relevant stability criteria for the model in terms of 
requirements on the parameter n: 

1. Non-negativity of the Gravitational Coupling: 

1-- 72n > 0, => n< — . 

- 72 

2. Absence of Tachyon Instability: 

1 - 72n > 0, n > . 

3. Ricci Stability: 

n{12n- 11) > 0, :^ n> 



11 

12 



It is unfortunate for the quadratic theories that are power- 
law in the GB invariant, that the above requirements can not 
be jointly met, so that these theories are not theoretically con- 
sistent in the sense discussed above. For negative n < 0, in 
particular, the bounds 1-3 reveal that, while the effective gra- 
vitational coupling Geff can be a positive quantity, as long 
as, for negative n = —n (n > 0), 1 + 72n > always, the 
theory ([T]) develops a tachyon instability since the mass squa- 
red of the propagating scalar degree of freedom is a negative 
quantity (recall that we are considering just positive curvature 
backgrounds): 



72n 



216n 



^o<0 



As already discussed, this is a serious objection against the 
quadratic modification of gravity depicted by the action ([T]) 
with positive and negative n-s. For f{G) in ([T]) with inver- 
se powers of the GB invariant (negative n < 0), the model 
considered here is a particular (spin-2 ghost-free) case of the 
model studied in Ref. ITolfTlll . 

We want to comment that, for positive n > 0, the cons- 
traints established on dU by consistency with solar-system test 
and with present cosmic speedup lead to n < 0,074 11511 . This 
is an example of a toy model that passes the cosmological tests 
while being theoretically inconsistent. 

Our previous results also apply to the more general model 
f{G) = mI^I" of Ref. H, since, in this case. 



' This is true only for power-law functions of the GB invariant. 



h=ntiG\Gr\ kG = 

so that, at constant curvature. Go 
positive magnitude. 



n{n-l)^AG^-\ 
= \G\ = X'^/&, is always a 



IV. INSTABILITIES IN THE TOY MODELS OF REE (H 

The toy models ^j^ were proposed in Ref. lUTIl on the ba- 
sis of several mathematical requirements on the function f{G) 
that make them cosmologically viable models. The most cru- 
cial condition to be satisfied is fgg > 0, which is required to 
ensure the stability of a late-time de-Sitter solution, as well as 
the existence of standard radiation/matter dominated epochs. 
The explicit f{G) models (|2]|4l) show a stage of cosmic accele- 
ration followed by the matter era. Although these models have 
been shown to be free of ghost instabilities, on the light of the 
discussion in the former sections we shall explore them with 
regard to the stability of vacuum solutions of constant curva- 
ture, by focusing on the requirements of Ricci stability and of 
non-negativity of the effective gravitational coupling. 



A. Toy model given by equation (|2) 



The first check will be to ask for positivity of the effective 
gravitational coupling SnGeff 



fi^a , 1. e.. 



1 



4X3/0^ > 



Let us write the function f{G) given by equation ^, in the 
following form: 

/(e) = AV^Karctane - ^ ln(l + e) - «] , 

where we have introduced a new field variable ^ = Q//3 {^q = 
X2/6/3), so that /e = (3fg, and /^^ = /J^/ee- Hence, 

The requirement of non-negativity of Ge// becomes into the 



fe = AV^arctan^, /^^ 

;ment of m 
following inequality: 

X^ - 144A V^X^ + 36^2 > . 

Since A (also /3, since the square root y/]3 has to be real) is 
non-negative, hence there is an interval in the values of the 
curvature Xq lying between the smaller and larger positive 
roots of the polynomial equation Xq - IMXy/^X^ + 36/3^ = 
0, where the above bound is not satisfied. Here a cautionary 
note is necessary: the value of the background constant curva- 
ture Xq can no be chosen at will, it is completely determined 
in terms of the free constant parameters A, /?, a, by the tra- 
ce equation (fT9] l which, in the present case, amounts to the 
following non-algebraic equation: 



In 1 



11 
36/32 



AV^ 



2a 



Since the LHS of this equation is always positive, then, any 
root Xq of this equation will obey Xq > 2a\y/]3. 

That the model does not develop gravitational ghosts is evi- 
dent from the following inequality, expressing the requirement 
of Ricci stabihty (/^ + 2Xlf^g > 0): 



arctan 



72/3^2 



6/3 y ^ 36/32 + X* 



>0 



which is always met. The model is also spin-0 tachyon-free, 
since the constraint fgg > was one of the requirements 
imposed in lUTIl to obtain the viable models (|2l|4li. 



B. Toy model given by equation (|3) 
In this case, since 

f{g) = —=Q arctan f -^ j - a\^ , 

or, in terms of the variable ^ = Q//3 

fiO = A V^(^ arctan f~ a) , 
then 



/c = Avp I arctan^ H — 



? 



e 



/ee 



_2WP_ 



The trace equation ( fT9l ) leads to the following 5th-order alge- 
braic equation to determine the background curvature Xq as 
function of the free parameters a. A, /3, 



X^ 



Xo 



36/32 + x^ 2\y/p 

Non-negativity of the effective gravitational coupling can be 
expressed in the form of the following inequality: 

X^-ZQ- UZP^Xq + 8 • m'^aXp^''^ > , 

where the trace equation has been taken into account. Extrac- 
ting useful qualitative information from the above inequality 
is very difficult. Absence of Ricci instability (fg + 2Xg /g^ > 
0) is guaranteed since (/3 > 0): 



arctan — 2- + — ^^- ,.°LJ:J > 



X'^ 



6/3 



6/3(1 +^4/36/32)2 



C. Toy model given by equation (|4) 
According to (|4| we have: 



.f{g) = X^l31n 



cosh 



/3 



»Av/^ 



For simplicity of writing let us to introduce the following pa- 
rameters: 



t = ^ = M Wo 
^ " /3 6/3 ' ° 



sinh^ 
cosh ^ 



then, non-negativity of the effective gravitational coupling 
Geff and Ricci stability requirement, can be written in the 
form of the following inequalities: 

^o'(l - W^o ) < ^, PWo + 2Xo2(l - W^) > , 



respectively. While the latter requirement is always fulfilled, 
non-negativity of Geff requires a careful examination in pa- 
rameter space due to the bound X§{1~ W§) < /3^/^/4A (the 
trace equation has to be considered). 

Although the models (|2l[3]i of reference lUTIl seem to be 
theoretically consistent also from the point of view of stabi- 
lity of vacuum, the cosmological constant term may be nee- 
ded in general for the cosmological viability of these models 
(a 7^ 0). As the authors of lUTIl recognize, the need of such 
a constant term makes the models (|2][3]) less attractive from a 
theoretical point of view. 



V. INSTABILITIES IN THE TOY MODELS OF REF. ^ 

Models (|5]|7]i where explored from the point of view of their 
asymptotic properties in Ref. JlSll . The authors demonstra- 
ted that cosmologically viable trajectories in the phase spa- 
ce of the model can be found. In particular a stable de Sit- 
ter point, corresponding to late-time cosmological evolution 
- co-existing together with phantom like stage -, is found. 
The behavior of the relevant cosmological parameters of ob- 
servational interest is quite similar to the one obtained in the 
ACDM-model. All of these models produce an elongated mat- 
ter dominated epoch followed by a radiation dominated epoch 
lUal . On the basis of their study the authors concluded that 
the toy models (|5]|6| are cosmologically viable. Here we aim 
at exploring the above toy models from the perspective of the 
vacuum structure of their linearizations in the sense discussed 
in the former sections. 



A. Model given by equation (|5) 
This toy model is given by the function: 

f{g) = ag'/^ + pg'^\ 

where a and (3 are arbitrary constants. In this case the trace 
equation ( fT9] l is easily solved to yield 



63/4/3 



4(1 + a/ V6) 



Hence, for instance 



i/2_4(V6 + a)2 i/4__2(\/6 + a) 



9/32 



3/3 



etc. Besides 



fg^^a + V6f{a~Vy2), 



and 



JQQ 35/36 



{a + V&f{a^^2ji) 



Non-negativity of the effective gravitational constant Geff'- 
1 — 'i^ofgg > 0, is achieved whenever 

17 [2 

while, absence of tachyon instability {Geff > 0, fgg > 0), 
imposes the following constraint; 

-{a + V6f{a - 7273) > ^ i^y^ < a < y^ . 
Ricci stabiUty: /^ + 2X1 j^g > 0' requires that 



33/32 



[oi + V6)2(8V6 - 25a) > ^ a < 



24 
25 



We see that Ricci stability, non-negativity of the effective gra- 
vitational coupling, and absence of tachyon instability, can be 
simultaneously fulfilled only if 

17/2 24/2" 

18V3-"-25V3' 

or if: 0,7711 < a < 0,7838. By comparing this resuh with 
the best-fit value of a reported in Tab. I of reference Ill9ll : 



a = 0,00084 



+0,00016 
-0,01632 ' 



one sees that this value (including error bars) does not fit in- 
to the above narrow interval where the model is theoretically 
consistent. 

The conclusions drawn from the above results suggest that 
cosmological viabihty of the model (|5]) is in serious trouble if 
one requires also its theoretical consistency. 



B. Model given by equation (|6) 

Another of the toy models proposed in lUSll . whose cosmo- 
logical viability was demonstrated in Ref. 11911 due to its good 
fitting to combined cosmological data sets, is given by the fo- 
llowing function: 



Hence, 



Zee 



j{Q) = a{g^i' - /3)2/3 



k = \{Q'^'-fir'i'Q-^^\ 



1(^3/4 _/3)-4/3 (^3/4 _^/2)g-5/4 



The trace equation for this model (considering, as we have 
done, a vacuum background of constant curvature where Qq = 
Xq/6) can be written as: 



where we have defined ^ = P/Qq = (S^^'^P/ Xj" . The only 
real root of the latter equation is given by 



i-l^rl, r.^[L 



1/3 



a 



{61)^3 



with ^ = 9 + a/81 + 36V6/a^ (since I is real, then, a < 
— 1,02872). In terms of the variable ^, taking into account the 
trace equation above, one has that: 

■'^ 2 \I3J 2C-1 ' 



and 



Jag — 



V6 ei2-o 



'^^ 8/32(2e-i)(e-i) 

Hence, non-negativity of Geff will imply that 



3 _ 2e + 15g - 35 



>0 



which leads to 

-oo < C < -9,3681 ; 0,5 < ^ < 1 , 
and, also 

1,8681 < ^ < oo . 
Absence of tachyon instability can be written as 



/< 



_V6 e{2-o 

^^ 8/32(2e-l)(e-l) 



>0 



e<l/2, l<e<2 



Additionally, Ricci stability bound is fulfilled if 



2/3 



5-2C 



/?; (2e-i)(e-i) 



>0 ^ e< 1/2, 1<C< 5/2 



Combining the above stability requirements for theoretical 
consistency of the model ^, yields to the following cons- 
traint on the parameter ^, that has to be fulfilled by the model 
if it is to be considered as theoretically consistent: 

^ >0; 1,8681 <C<2, 

while, for negative /3-s: 

/3 < ; -oo < ^ < -9,3681 . 

If one substitutes these limits back into the trace equation, one 
gets the corresponding allowed ranges of the parameter a: 



/3 > ; -1,6984 <a< -1,6706 



13 <0; -1,9764 <a< -1,02872 



and 



On the other hand, in Tab. I of Ref. 111911 the reported best-fit 
value of the parameter a: 

,+0,000083 



-0,00014" 



0,001439 ' 



SO that, considering error bars, the model is only marginally 
consistent. 



C. Model given by equation 

The last toy model we will explore here is the one given by 
the following function: 



f{G) = aVOexp ( - 



It will be useful to introduce the parameter ^ = /3/Qq, which, 
depending on the sign of /3 can be either positive or negative. 
We have 



/o = "Wt (^ ' /( 



'-''-Wl 



^^1-20 



and 



In terms of these parameters the trace equation ( [19] ) can be 
written as 



V6 



(1 + 20- 



This equation will be used to remove the exponential from the 
expressions below. 

Non-negativity of Geff can be written as the following 
constraint 



4^2 I 73 f _ 35 
'' ^ 18^ 36 



>o, 



1 + 2^ 

which can be solved to yield: 

-1,2141 < C < -0,5 ; 0,2002 < C < oo . 

Once Geff is checked to be non-negative, absence of spin-0 
tachyon mode amounts to fgg > 0: 



1 - 4C - 4^^ 

i + 2e 



>0 



Hence: 



-oo < C < -1,2071 ; -0,5 < ^ < 0,2071 



Joint fulfillment of the above requirements, meaning non- 
negative Geff and absence of scalar tachyon, is possible only 
if the parameter ^ falls within the following, very narrow in- 
tervals: 

-1,2141 < ^ < -1,2071 ; 0,2002 < ^ < 0,2071 . 

Ricci stability 



f° + 2X'f" >0 => 5 -^^^-^4^' 
Jg + ■^^o.lgg > U ^ T+^f^ 

imposes additional constraints: 



>0, 



-oo < C < -1,1052 ; -0,5 < C < 0,1885 
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The latter ranges match with the ones where joint fulfillment 
of Geff > and of f^g > 0, is allowed, only if -1,2141 < 
^ < —1,2071, or, in terms of the original parameters: 

-4,9706/3 < X^ < -4,9419/3 , 

which forces the parameter /3 to be negative. Besides, if one 
substitutes the values of ^ within the above range, back into 
the trace equation, one gets that the parameter a takes values 
within the following interval 

5,7752 <a< 5,7916 , 

which, including error bars, is inconsistent with the best-fit 
a-value 

a = -0,00012_o_oo295 > 

reported in Tab. I of reference lll9r . Therefore, observational 
viability of this model is obviously in conflict with its theore- 
tical consistency. 



VI. DISCUSSION 

A careful study of the stability of higher-order modifica- 
tions of gravity has been proved to be useful in ruling out 
theoretically inconsistent models. This is particularly true for 
/(X)-theories. It has been demonstrated, in particular, that 
these theories are equivalent to scalar-tensor gravity, with a 
curvature-dependent mass squared of the scalar field (see the 
review ll22ll and references therein). The dependence on the 
curvature makes it possible to explain the modifications in the 
large scale caused by the scalar mode: very light mass of the 
scalar implies that the modification of gravity has implications 
for cosmology, while, heavy mass in the neighborhood of ga- 
laxies and star systems, makes the scalar mode to be shielded 
at small scales thus rendering Newton's law of gravity valid. 
When dealing with arbitrary quadratic modifications of gra- 
vity of the form F{X, Y, Z), the situation is much more com- 
plicated. In this case the original theory is shown to be equiva- 
lent to a multi(quadruple)-scalar-tensor theory with the scalars 
coupled to the different curvature invariants [12]. This makes 
the equivalence useless to study the field content of the theory 
as well as to study observational constraints on it. In order 
to study such fundamental issues one, usually, limits oneself 
to considering excitations around a vacuum state (a maxima- 
lly symmetric spacetime with a constant curvature) (al- First 
one expands the action around spaces of constant curvature, 
and then one incorporates linear perturbations of the metric 
around vacuum solutions. The spectrum of perturbations is 
then revealed. For general theories F{X, Y, Z) there arise 8th 
degrees of freedom: two associated with the spin-2 (massless) 
graviton, one with the spin-0 scalar mode, and five with the 
spin-2 massive Weyl state JtI-IqI [l2ll . The latter happens to be 
a ghost due to the negative sign of the Weyl term after the 
linearization. Much information about the original (unpertur- 
bed) theory can be extracted from the properties of the above 
spectrum of perturbations around vacuum background spaces. 



Actually, if one wants to have a consistent expansion of the 
theory around such vacuum states, then one has to require, in 
the linearized theory: i) the effective gravitational coupling to 
be non-negative, so that the massless graviton is a normaliza- 
ble state, ii) the mass squared of the scalar propagating mode 
to be non-negative since, otherwise, the spin-0 mode would 
be a tachyon state leading to disastrous destabilization of the 
vacuum, and iii) the negative mass squared of the spin-2 Weyl 
mode to be large enough as to decouple from the linearized 
spectrum (no unphysical Weyl ghost, or "poltergeist"). The 
last requirement (spin-2 ghost free theory) is achieved if one 
considers particular quadratic modifications of the kind stu- 
died here: F{X, g) ^ X + f{Q). At the Hnearized level these 
are scalar-tensor theories of gravity. 

Following discussions within /(,Y)-theories, one might in- 
voke other instabilities that are of importance in cosmological 
applications. In particular, the so called Ricci, or also, Dolgov- 
Kawasaki instability, associated with the fact that the graviton 
might become a ghost if the curvature increases |3, 22-24fl. 
As it is usually done in cosmological applications of f{X)- 
theories, working with F{X, Y, Z) models one can write the 
LHS of the field equations in the form of standard Einstein's 
equation, with a non-standard RHS. In this case, the effective 
gravitational (Newton's) constant G^ depends on the curva- 
ture - strictly speaking, it depends on the X-derivative of the 
function F{X, Y, Z): G"^^ (x F^^ -, so that it varies with the 
course of the cosmic expansion. It might happen that, as the 
cosmic expansion proceeds the value of G^ increases. This 
would has a disastrous impact on the theory, since the curva- 
ture might grow without limits 11221 uM . Therefore, one has to 
require the X-derivative of the effective Newton's constant to 



be non-positive: dG^^ /dX < 0, which means that Gj^ is 
a non-increasing function of the background curvature. This 
will warrant that the graviton does not become a ghost as the 
background curvature grows. 

A combination of the above requirements (non-negativity 
of the effective gravitational coupling, Ricci stability, absen- 
ce of scalar tachyon) might be very useful to rule out ap- 
parently viable quadratic theories of gravity. In this paper 
we have checked several such cosmologically viable theo- 
ries iH [H [l7l4ll, where the function F{X,Y,Z) = 
F{X, Q) = X + f{G), regarding the above discussed sta- 
bility tests. It arose that some of these "viable" models did not 
pass the stability tests considered here. In particular, if one ta- 
kes into account the best-fit values that make possible to talk 
about cosmological viability of some of the models, one sees 
that, either the models clearly fail to meet the stability requi- 
rements, or they are on the edge of (not) being theoretically 
consistent. 

It has been demonstrated that theories based on Lagran- 
gians which contain powers of the GB invariant (for negative 
powers n < 0, these are a particular case of the model pro- 
posed in Ref. lllOll ). are Ricci unstable and, also, develop sca- 
lar tachyons. Ricci and tachyon instabilities combined, provo- 
ke gross destabilizing effects with catastrophic consequences 
for the theory. This result applies also to models of the kind 
f{g) = /i|g|", thoroughly explored in Ref. 03 . A careful 
study of the models (|5]|7]i [iJI], reveals that, when taking in- 
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to account the reported best-fit values of the parameters of the 
models, the stability criteria are either clearly not met (models 
(|5]l and (Q), or are at the edge of being not met (model ^). 
While the situation with models (|2l|4|l jlTT is better, neverthe- 
less, careful consideration of the free parameters of the theory 
is required. Besides, even when these models were theoretica- 
lly consistent, the cosmological constant term may be needed 
in general for their cosmological viability {a 7^ 0). As the aut- 
hors of lUTIl recognize, the need of such a constant term makes 
these models unappealing from a theoretical point of view. 



VII. CONCLUSIONS 

In this paper we have performed a careful study of the stabi- 
lity of vacuum linearization of several cosmologically viable 
models based on quadratic modifications of gravity. This kind 
of study is complementary to other existing tests such as the 
application of the dynamical systems tools Iil8i1 , and model- 
fitting of the combined data sets II19I1 . We have demonstrated 
that several models that pass the latter tests, fail to have a con- 
sistent vacuum linearization due to the instabilities caused by 
the presence of ghost and tachyon states. On the basis of the 
mere occurrence of any of these instabilities a given, suppo- 
sedly viable cosmological model, might be rule out. 

We have paid special attention to: i) the development of 
ghost graviton, associated with increasing values of the effec- 
tive Newton's constant G^ with the background curvature 
(the so called Dolgov-Kawasaki instability, also called Ricci 



instability), ii) to the requirement of non-negativity of the ef- 
fective gravitational coupling Geff, ensuring that the graviton 
is not a ghost state, and iii) to absence of a destabilizing spin-0 
tachyon perturbation. 



The results of our study suggest that models ^j/Q} llisll ex- 
hibit vacuum instabilities of the kind considered if, at the sa- 
me time, the best-fit values of the free parameters taken from 
III9I1 are considered. While models ^j^ are less problema- 
tic in this respect - although a careful study of the parameter 
space is necessary -, the models containing powers of the GB 
invariant lllC)l[l4ll are clearly to be rule out. 

Our conclusion is that the inclusion of quadratic modifica- 
tions of gravity, the way these modifications have been consi- 
dered so far, results in very unappealing models. Perhaps an 
alternative simple way to introduce these modifications is th- 
rough the Dirac-Born-lnfeld deformation strategy, according 
to which the Lagrangian density C = \/\g\L is to be replaced 



by £ = ^y\g\e\{^l + 2eL/\ - 1) (e = ±1, A is an addi- 
tional energy scale). This strategy has been invoked to erase 
unpleasant singularities from the original theory. In the con- 
text of the quadratic theories of gravity it might help to avoid 
or smooth out several of the instabilities considered here luOll . 
This will be the subject of future publications. 
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